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We present exemplary results of extensive studies of mechanical, electronic and transport properties of covalent 
functionalization of graphene monolayers (GML) with -NH 2 . We report new results of ah initio studies of 
covalent functionalization of GML with -NH2 groups up to 12.5% concentration. Our studies are performed 
in the framework of the density functional theory (DFT) and non-equilibrium Green's function (NEGF). 
We discuss the stability (adsorption energy), elastic moduli, electronic structure, band gaps, and effective 
electron masses as a function of the density of the adsorbed molecules. We also show the conductance and 
I(V) characteristic of these systems. Generally, the stability of the functionalized graphene layers decreases 
with the growing concentration of attachments and we determine the critical density of the molecules that 
can be chemisorbed on the surface of GLs. Because of local deformations of GLs and sp 3 rehybridization of 
the bonds induced by fragments, elastic moduli decrease with increasing number of groups. Simultaneously, 
we observe that the functionalizing molecules stretch the graphenes lattice, the effect being more pronounced 
for higher concentration of adsorbed molecules. We find out that the GLs functionalization leads in many 
cases to the opening of the graphene band gap (up to 0.5302 eV for 12.5% concentration) and can be therefore 
utilized in graphene devices. The new HOMO and LUMO originate mostly from the impurity bands induced 
by the functionalization and they exhibit parabolic dispersion with electron effective masses comparable to 
ones in silicon or gallium nitride. 
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I. INTRODUCTION 

Nowadays graphene, a 2-dimensional semimetalic 
monolayer, a form of sp 2 hybridizated carbon, attracts a 
lot of research activity owing to its unique electronic, me- 
chanical and thermal properties- 1 . Graphene layers (GLs) 
are emerging as the very promising candidates for a new 
generation of electronic devices 2 - - — . However, the GML 
has zero energy band gap. This hinders direct applica- 
tion of graphene layer in field effect transistors (FETs) 
and the functionalization of GLs could be a remedy for 
this problem. The possibility to generate controllable 
band gap in graphene without significant deterioration 
of the remaining advantageous properties is desirable. 
Furthermore, covalent functionalization enables usage of 
graphene monolayers as sensors^. 

We have performed extensive studies of the stability, 
electronic and transport properties of the functionalized 
graphene with various molecular groups^. In this article 
we focus on the problem how the properties of GMLs 
change after the functionalization with -NH2 . 

II. CALCULATION DETAILS 

Our studies of GML functionalized with -NH2 are 
based on ab initio calculations within the framework of 
the spin polarized density functional theory (DFT)2ii2 
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as implemented in the SIESTA package 11,12 . The so- 
called PBE form of the generalized gradient approxima- 
tion (GGA)^ 3 - is chosen for the exchange correlation den- 
sity functional. The computations involving full geom- 
etry optimization have been performed employing fol- 
lowing parameters of the SIESTA package that deter- 
mine numerical accuracy of the results: double-£-plus- 
polarization basis, kinetic energy mesh cutoff of 350 Ry, 
the self-consistency mixing rate of 0.1, the convergence 
criterion for the density matrix of 10~ 4 , maximum force 
tolerance equal to 0.01 eV/ A. Calculations were per- 
formed for supercell geometry with graphene layers sep- 
arated by a distance large enough to eliminate all kinds of 
spurious interactions. For the band structure and density 
of states calculations, the Mesh Cutoff has been increased 
up to 550 Ry. For better imaging of density of states, the 
peak width for broadening the energy eigenvalues has 
been set to 0.05 eV. To investigate the various densities 
of the functionalizing amines, we have taken large lat- 
eral unit cells that contain 4, 9, 16, 25, and 36 standard 
two-atomic graphene unit cells arranged respectively in 
the patterns (2x2), (3x3), (4x4), (5x5), (6x6) (i.e., these 
supercells contain 8, 18, 32, 50 and 72 carbon atoms, 
respectively) and placed a functionalizing group in the 
enlarged unit cells. 

The adsorption energy, which has been calculated ac- 
cording to the formula: 

E ads = (Etot (MGL + NH 2 ) - E tot {MGL) - E tot (NH 2 )) 

(1) 

is used as the measure of the stability of the systems 
studied. E to t{MGL + NH 2 ) is the total energy of func- 
tionalized graphene, E to t(MGL) is the total energy of 
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graphene monolayer, and E to t(N H2) is the total energy 
of the amine group. 

Elastic properties such as Poisson ratio, Young's, shear 
and Bulk moduli have been also calculated. The Poisson 
ratio values, has been obtained according to the equation: 
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where Aa(b) is change of lattice constants in strained 
GML. The most interesting quantity, Young's modulus, 
have been determined - from components (era) of the 
stress tensor: 
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where en is strain. Bulk modulus (Eq. 0]) and Shear mod- 
ulus (Eq. [S|) , which could be easily derived from (Eq. [2]) 
and (Eq. [3]) are given by the following mathematical for- 
mula: 
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We have used TranSIESTA 14 to study electronic co- 
herent transport in the systems employing NEGF tech- 
nique, within the Keldysh formalism, a rather standard 
procedure for treatment of coherent transport 15 . The 
structures have been treated as two-probe systems with 
the central scattering region sandwiched between semi- 
infinite source (left) and drain (right) electrode regions 
(as depicted in Fig. EJa)). The conductance can be ex- 
pressed as follows^ 5 -: 



G(E) = G Tr[T R G r T L G% 



(6) 



where G — 2^- is the unit of quantum conductance 
and G r ( a ) is retarded (advanced) Green's function. The 
current through the scattering region has been calculated 
according to Landauer-Buttiker formula, i.e., assuming 
the limit of small bias, which is justified for the range of 
the external voltage considered in the present study: 
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where Vds is equal to the electrochemical potential dif- 
ference between the left and right electrodes (eVcs^ML — 
Hr) and T is the transmission. 



III. RESULTS AND DISCUSSION 

Amine groups bind covalently to the graphene layer, 
implying sp 2 to sp 3 rehybridization and local deforma- 
tion of the graphene plane, which are well known ef- 
fects observed in chemical doping graphene layers and 




FIG. 1. (a) The optimized graphene monolayer functionalized 
with one -NH2 group per (3x3) supercell. (b) Simulated STM 
image of structure visualised above (Vj,i as = lV). Also super- 
posed is a ball-and-stick model of the functionalized graphene 
lattice. 



nanoribbons 5 -^ - — . C-C bond length in the closest prox- 
imity of the functionalizing group increases to 1.51 A, 
whereas C-C bond lengths away from the group are very 
close to the bond length in pure graphene (1.43 A). The 
carbon atom that is directly bonded to the group sticks 
out from the graphene sheet (see Fig. Q] (a)), the effect 
being more pronounced for higher concentration of ad- 
sorbed molecules. The STM image depicted in Fig. LTJ 
(b) is calculated under the Tersoff-Hamann theory^ and 
visualized using WSxM 20 . It shows bright spots corre- 
sponding to amine groups above graphene lattice. It is 
clearly seen that amines break the hexagonnal symmetry 
of graphene layer and create scattering centers. This is 
consistent with recent experimental worki£. 

The adsorption energy, used as a measure of the stabil- 
ity, decreases with growing concentration of groups. Be- 
cause of local deformations of GLs and sp 3 rehybridiza- 
tion of the bonds induced by fragments, elastic moduli 
diminish with increasing number of groups. The high- 
est reduction is equal to 15.5% in Shear, 14% in Young's 
and 9.5% in Bulk modulus (see Tab. [I]). In comparison to 
standard materials, the elastic properties of amine func- 
tionalized GMLs are still superior and its potential usage 
in solar cell is secured 5 -. 

We find out that the GLs functionalization leads in 
many cases to the opening of the graphene band gap (up 
to 0.5302 eV for 12.5% concentration - see Tab. HJ and 
can be utilized in graphene devices. 7 Similar results, but 
for bilayer graphene, were obtained by Boukhvalov^. 

The fundamental band gaps for the functionalized 
structures generally inrease with increasing concentra- 
tion of amines and typically are indirect. It is specified 
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graphene+NH 2 




FIG. 2. Band structure and DOS of graphene functionalized 
with amine molecule (one per (3x3) supercell) are plotted with 
black lines, whereas these quantities for pristine graphene are 
plotted with blue dotted lines. Inset: zoom of band structure 
around Fermi level (shifted to 0) along direction M— >K in the 
Brillouin Zone. Indirect band gap is marked by the yellow 
dotted arrow. 

in Fig. [H The new HOMO and LUMO originate from 
the impurity bands induced by the functionalization and 
exhibit parabolic dispersion. Those impurity levels re- 
sult mostly from nitrogen atom and 7r-bonds. The effec- 
tive electron masses (m e ff) and mobility depend rather 
strongly on the concentration of dopants. 

At some concentrations, they are comparable to masses 
in silicon or gallium nitride (see Tab. H]). For (3x3) and 
(6x6) supercells m e // are close to 0, which is consistent 
with the phenomenon of zero band gap for supercells 
with dimensions being multiples of three graphene lat- 
tice constants described by Garca-Lastra^. Breaking of 
the hexagonal symmetry of the layer is also clearly seen in 
effective masses. They are typically slightly bigger along 
K-r than K-M direction. 

From the transport calculations (see Fig. [3]) , one can 
see that this type of functionalization hardly changes 
the conductance (Fig. EJe)) and I(V) characteristic 
(Fig. Eld)) in comparison to pure graphene monolayer. 
Nevertheless, it is clearly seen that the conduction 
around Fermi level decreases with increase of concentra- 
tion of adsorbants. The critical voltage for the non-zero 
current is the lowest for the system with two amines and 
the highest for the pristine graphene. Also for low bias, 
the current for functionalized systems is larger than for 
unfunctionalized one. However, for the bias larger than 
0.6 V, we observe the largest current for pristine graphene 
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FIG. 3. Three systems studied for transport properties: pure 
graphene monolayer (a), graphene functionlaized with one (b) 
and two (c) NEb groups. Electrode regions contain 16 carbon 
atoms, whereas scattering regions contains 80 carbon atoms 
and zero, one, and two -NH2 groups for cases (a), (b), and (c), 
respectively. Periodic boundary conditions are conserved in 
perpendicular direction to the z axis. The voltage is applied 
along the z axis, (d) I(V) characteristic for pure (blue dotted 
line) and functionalized graphene with one (black solid line) 
and two (red solid line) -NH2 groups, (e) The transmission 
spectra for three different source-drain voltages calculated for 
three systems described above. 

and the smallest for the system with higher concentration 
of functionalizing amines. Let us remark that the conduc- 
tance of graphene functionalized with amines drops much 
slowlier with concentartion of dopants than in the case of 
functionalization with hydroxyl or carboxyl groups. It is 
mostly due to the fact that nitrogen acts as donor. This 
is consistent with the experimental results of Barakel^. 
The increasing concentration of primary amines cova- 
lently bound to graphene monolayer causes increase in 
the chemical reactivity of the surface, while the electri- 
cal conductivity decreased. However, even highly ami- 
nated graphene, up to 20%, is conductive enough to be 
be used for DNA detection as bio-attachment platform 
in a biologically active field-effect transistors. Therefore, 
it is plausible that these functionalized structures can be 
commercially used as biosensors 5 - . 



IV. CONCLUSIONS 

In summary, we have performed ab initio calcula- 
tions of the mechanical, electronic and transport prop- 
erties of the graphene monolayer functionalized with - 
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TABLE I. Graphene monolayers functionalized with amines: c- concentration of groups, E ads - adsorption energy, Y- Young's 
modulus, S - Shear modulus, B- Bulk modulus, r\ - Poisson ratio, E gap -fundamental band gap, m^ M and - effective 

mass of electron in K point along the direction K-M and K-P, respectively. 



supercell 


c[%] 


E ads [eV] 


Y [TPa] 


S [TPa] 


B [TPa] 




K-M 
m e// 


[m e ] mf,./ [m e ] 


(6x6) 


1.39 


-0.883 


1.02 


0.46 


0.43 


0.11 0.0058 


0.005 


0.010 


(5x5) 


2.00 


-0.846 


1.03 


0.47 


0.44 


0.11 0.0732 


2.675 


2.774 


(4x4) 


4.35 


-0.825 


1.01 


0.45 


0.43 


0.11 0.0852 


3.368 


6.144 


(3x3) 


5.56 


-0.799 


0.98 


0.44 


0.42 


0.11 0.1535 


0.038 


0.039 


(2x2) 


12.5 


-0.739 


0.90 


0.40 


0.42 


0.14 0.5302 


20.148 


20.720 



NH2 molecules. Covalent functionalization changes the 
local and global structure of the functionalized system, 
causes sp 2 — > sp 3 rehybridization, reduces the elastic 
moduli, and introduces impurity levels around Fermi en- 
ergy, which allows us for band gap engineering. This 
type of functionalization slightly reduces current and 
only moderately modifies the I(V) characteristics in com- 
parison to the pristine graphene. 
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